Phase diagram of imbalanced strongly interacting fermions on a one dimensional 

optical lattice 
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We show that the Hubbard Hamiltonian with particle-assisted tunneUng rates -recently proposed 
to model a fermionic mixture near a broad Feshbach resonance- displays a ground state phase 
diagram with superfluid, insulating, and phase separated regimes. In the latter case, when the 
populations are balanced the two phases coexist in microscopic antiferromagnetic domains. Macro- 
scopic phase segregation into a high-density superfluid of molecules, and a low-density Fermi liquid 
of single atoms appears in the density profile above a critical polarization pc- 
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I. INTRODUCTION 

The recent experimental observation of normal- 
superfluid transition in two-component Fermi gases [J, [H, 
[3] has attracted a great deal of attention on the subject 
of strongly interacting confined fermionic atoms. In this 
context, an important result consisted in the observation 
on the density profile of phase separation (PS) between a 
superfluid core and an external polarized shell, depend- 
ing on the population imbalance p and the interaction 
strength across a Feshbach resonance [1, [5|. The theo- 
retical investigation of PS in the framework of uniform 
Fermi gases [6] provided an efficient description of the 
phase diagram along the BCS-BEC crossover. In partic- 
ular, homogeneous phases of fully polarized and partly 
polarized normal gases or of unpolarized and partly po- 
larized superfluids have been recognized, as well as the 
coexistence of some of these phases by varying the gas 
polarization and the interaction strength between atoms. 

A description of the above phenomena when the 
fermionic atoms are confined by anisotropic optical lat- 
tices would represent a major achievement. On the the- 
oretical side, it has been speculated that, with respect 
to real materials, these systems could provide a much 
neater realization of the Hubbard model 0, @] • The very 
recent experimental observation [9*, of the formation 
of an insulating phase in a repulsively interacting two- 
component Fermi gas on an optical lattice confirms this 
hypothesis, thus, enabling the use of such systems as a 
laboratory to investigate the onset of many phenomena 
in real materials, one for all high temperature supercon- 
ductivity. 

There is still no full agreement on the fact that 
Hubbard-like model Hamiltonians can support evidence 
of macroscopic phase segregation. Moreover, it has been 



observed that the Feshbach resonance, besides inducing 
the desired strong interaction among atom s [ll| . can 
cause highly nontrivial effects [l2|, [3, Il5|. which 
should be effectively included in the lattice model. For 
instance, a broad resonance -typical of the experimental 
setups with fermionic atoms- could generate a multiband 
distribution of the atomic population, which in turn af- 
fects the effective correlation between atoms on neigh- 
boring sites. In such a case, it was shown [15] that a lat- 
tice resonance model is the generic model which would 
proper ly describe the system. The latter can be mapped 
0, i^^^ extended Hubbard model with correlated 
hopping, previously known in literature as the Simon- 
Aligia model Hamiltonian [18]. 

In this work, we provide exact analytical and numerical 
evidence that in the regime of strong interaction between 
fermions the Simon- Aligia model reproduces the scenario 
expected for ultracold fermionic gases, exhibiting normal, 
superfluid, insulating, and PS regimes. Based on such re- 
sults, we also propose a simple explanation of the micro- 
scopic mechanism driving the transition to macroscopic 
phase segregation: while phase coexistence is present for 
appropriate parameters in a large range of filling values, 
the short range antiferromagnetic order masks this fea- 
ture when the two populations are balanced and only for 
a critical polarization pc macroscopic phase segregation 
appears. Above it, a regime of breached pairs is entered. 



II. LATTICE MODEL HAMILTONIAN 

In the case of fermionic atoms with two internal states 
on a one-dimensional optical lattice, the Simon-Aligia 
model can be written as [16,] 



HsA = 



E 



[t + 5g{ni-cr + ^j-cr) + Stni-fjTij- 



(1) 
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where a = {t, i} identifies the two internal states of the 
fermionic atoms {—a standing for the opposite of a), (z, j) 
denotes two neighboring sites, and rii = n^^ + riii. Here 
t describes the direct hopping of atoms of a given pop- 
ulation between neighboring sites, while Sg = g — t and 
St = t -\- tad — "^g are the deviation from the direct hop- 
ping case (in which Sg = = St)^ induced by correlations 
in proximity of a wide Feshbach resonance. More pre- 
cisely, g describes the configuration tunneling, in which 
one atom is transferred onto a site already occupied by 
an atom of different specie: it amounts to coupling two 
fermionic atoms into a dressed molecule, tad accounts 
for the motion of one atom between two already oc- 
cupied sites: de facto it exchanges a dressed molecule 
and a fermionic atom located on neighboring sites. Fi- 
nally, A is the energy cost of the dressed molecule, which 
works as an effective detuning parameter. Since in Eq. 
([1]) the number of atoms of both populations are con- 
served quantities, one can work at arbitrary fixed aver- 
age number of atoms per site (filling) n = N/L with 
= (Xli^i) (0 ^ ^ ^ 2), and population imbalance 
P={Einn-n,l)/N {0<\p\<l). 

In the context of highly-correlated fermionic materials, 
HsA has been widely studied [Tol, [2Q| , I2H . In particular, 
at p = it was recently shown [20] that -for the choice 
Sg = —t (i.e., g = 0)- its exact ground state phase dia- 
gram can be obtained, and it has been explored at n = 1. 
The structure of the latter is reminiscent for some aspect 
of the mentioned physics of cold fermionic gases, pro- 
viding evidence of both PS and an insulating behavior 
for appropriate values of interaction parameters. Both 
features were confirmed by numerical analysis also for 
< g < —tad [HI. Here we recast the exact solution of 
Ref. |20| to the present context, allowing for arbitrary p, 
and we explore the role of population imbalance on the 
density profile at and away from the exact case. 

III. RESULTS 
A. ^ = 

We first consider the ^ = case. A convenient rep- 
resentation of the model Hamiltonian ([1]) is obtained 
by introducing on-site Hubbard projection operators. 
These are defined as X^^ = \a)i{P\i, where \a)i are 
the states allowed at a given site i, and a = 0,t,|,2 
(|2) = I ID). In terms of these operators the Hamilto- 
nian H = HsA^Sg = —t) reads as (up to constant terms) 

H=Y. [t^f^T + tadX^^Xf] + A ^ Xf . (2) 

I 



egs{p, n) = min <--[(!- Ih) sin 



From a mathematical point of view, the choice Sg = —t 
implemented in Eq. ([2]) implies that also the total num- 
ber of dressed molecules Nd = (J^^-^P) is a conserved 
quantity. Hence, the ground state ojf H must be searched 
among the eigenstates of just its first two terms, namely, 
those with coefficients t and tad- These terms are de- 
generate with respect to the population of the single 
atoms and their eigenstates do not depend on p. One 
can show that in the thermodynamic limit, the ground 
state \^l^Gs{'^)) consists of a high-density core of length 
Lh with Nd paired atoms and Lh — Nd fermionic atoms, 
surrounded by a low-density shell of length L — Lh oc- 
cupied hy N — Nd — Lh single atoms, where L^ and Nd 
have to be determined variationally. Such a choice al- 
lows one to maximize the number of available low energy 
momenta. Explicitly, 



IV'Gs(n)) = \MN-Na-Lh))L-L,\MLh-Na))L, ■ (3) 



Here \^|Ja{Xf))L are the ground states of Nf effective 
spinless fermions on a L-sites chain, moving in a back- 
ground of empty sites a = / or dressed molecules a = h. 

For Lh = = Nd the ground state coincides with 
\i^i{Xs))L and describes a normal (possibly partly polar- 
ized) Fermi liquid of single fermionic atoms, which we de- 
note as Np. Moreover, for Lh = L^ the ground state be- 
comes \iph{L — Nd))L^ a (polarized) superfluid of dressed 
molecules and fermionic atoms {SFp). In particular, for 
Lh = Xd^ the latter consists just of dressed molecules, 
and we denote the state as unpolarized superfluid (SFq). 
Apart from these cases, \^pGs{'^)) is always characterized 
by PS. The two coexisting phases have different densi- 
ties n^, n/i, where can be expressed in terms of the 
variational parameters lh = Lh/L and rid = Nd/L as 



ni = l , n/, = 1 + — . (4) 

J- — lh 



The boundaries of the PS regions in dependence of 
the physical parameters can be derived by investigat- 
ing, at fixed p and n, the dependence on Ud and 
lh of the energy in the thermodynamic limit, namely, 
limL^oo ^{i)Gs{n)\H\il)Gs{n)). The ground state energy 
egs (p, n) is obtained upon minimizing the latter with re- 
spect to Ud and lh 



+ //,tadsin7r(2 - n/,)] + An^ ^ , (5) 
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where the expression in parenthesis is constrained by the 
actual value of the polarization since Ud < n/2{l — p). 
The optimal solution is obtained for p = 0, in corre- 
spondence with Ud = Ud and Ih = Ih- In order to dis- 
cuss the scenario of imbalanced populations, we notice 
that the minimization equation for Ud has formal solu- 



tion Ud 



ih 



arccos ■ 



2^a 



While at p = this sat- 



isfies the constraint for Ud^ in general it may miss such 
property for p 7^ 0. In fact, this happens above a critical 
polarization 

Pc = l-2fid/n , (6) 
when a regime of breached pairs is entered. The system 



At lower-enough attractive A value, the solution becomes 
fid < Ih- According to Eq. (|4j), this implies rih < 2: 
a regime is entered in which the superfiuid is polarized 
(SFp) since it contains 2 — Uh unpaired atoms and still 
coexists with Np. Such a region extends up to Ac, which, 
in the present case, is moderately repulsive. Since the size 
L — Lh of Np increases by enhancing A, also the value of 
the critical polarization pc above which the normal phase 
becomes fully polarized increases accordingly. Above it, 
the SFp phase reduces its size by breaking some pair: 
this case corresponds again to Ud < Ud- In both cases, the 
transition of the superfiuid phase in the PS regime from 
SFq to SFp is identified by the line Uh = 2. Moreover, 
the transition line to the breached pair regime (shaded 
region in Fig. [2|) is characterized by fid = n/2(l — p). 
Finally, the uniform phase of normal Fermi liquid (Np) 
is reached for A > Ac, the transition line being denoted 
by = T/i = 0. The physical parameters tad and n also 




Figure 1: (Color online) Number of pairs Ud versus polariza- 
tion p at A = 0, and: Sg = — t, St = 0.4t (continuous line, 
exact), Sg = — 0.8t, St = (red diamonds, numerical). Inset: 
critical polarization pc vs filling n. 



begins nucleating the normal phase within the superfiuid 
to accommodate the excess polarization and, correspond- 
ingly, the number of dressed molecules Ud diminishes with 
respect to the optimal value n^^, Ud = n/2{l — p) for 
P^Pci as shown in Fig. [TJ 

The ground-state phase diagram of the model de- 
scribed by as obtained from Eqs. (Hj)-® is reported 
in Fig. [2] in the A-p plane at a filling value n < 1 and 
t > ^ad > t/2. Six different regions can be identified, de- 
pending on the values of Ud and Ih in the solution. In the 
regime of deep attractive detuning at p = 0, the solution 
corresponds to the choice n^^ = = n/2: it describes 
an unpolarized superfiuid of pairs, SFq. In the same 
regime when p 7^ 0, the solution becomes Ih = rid < fid'- 
the ground state nucleates in S'Fq, a low-density normal 
phase of fully polarized Fermi liquid (Npp). In this case, 
PS is favored with respect to the polarized superfiuid, 
since tad ^ t; for tad > t the opposite scenario holds. 
When the detuning A is still attractive but moderate, as 
well as p, the solution reads fid = Ih < ri/2: the unpolar- 
ized superfiuid coexists with a normal Fermi liquid {Np). 



Figure 2: Exact phase diagram in the A-p plane at n = 0.9, 
Sg = — t, and St = 0.4t: the shaded region corresponds to 
macroscopic phase segregation in the density profile. 

play a relevant role in the above scenario. The filling n 
is crucial in the PS regimes, which is entered only for 
n > ni. In particular, the SFp + Np regime is the sta- 
ble ground states in the range ni < n < Uh^ whereas the 
SFo-\-Np state is stable up to n/^ = 2, ni and rih depend- 
ing on the physical parameters through Eq. (|4]). Also the 
uniform phase changes by varying n: while for n < 1 this 
is Np^ for n = 1 it becomes insulating, and for n > 1 it is 
SFp. A thorough investigation shows that by decreasing 
tad^ the phase diagram becomes more asymmetric with 
respect to the half filling configuration (not shown). 



B. gy^O 

We now discuss the g case. Since in this case no 
exact solution is available, one has to resort to numerical 
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simulations for finding the ground state of Hsa- At p = 
0, the situation was already investigated by means of 
the density- matrix renormalization-group algorithm [2l| : 
the essential features of the exact solution characterize a 
large range of g values, namely, < ^ < —tadi whereas 
a crossover region to the standard Hubbard regime is 
achieved when g > —tad and still tad < 0. The main 
difference with respect to the ^ = case is that the spin 
degrees of freedom become relevant. As a consequence, 
at p = when < g < —tad PS manifests itself at a 
microscopic level, by forming nanoscopic PS domains of 
size A (that turns out to be related to Ud) in which the 
coexisting phases rearrange their relative spins in order 
to maximize antiferromagnetic short-range order. 

Given the above interpretative framework, we thus ex- 
pect that p should not affect the results at ^ 7^ 
as far as the presence of phase coexistence is concerned. 
This is confirmed by the numerical data obtained by ex- 
act diagonalization on small clusters (L = 16) for rid and 
Pc dit g = 0.2t, compared in Fig. [T] with the exact re- 
sults at ^ = and same /7, t^d, and n values, showing 
accurate quantitative agreement. In both cases, increas- 
ing p amounts to force the orientation of an increasing 
number of single atoms, while Ud remains constant and 
equal to the p = optimal value up to p = Pc- Here, 
all the single atoms belong to one of the two internal 
states and, for p > Pc^ some pair must break with re- 
spect to the optimal number fid. What is expected to 
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Figure 3: (Color online) Static charge structure factor N{q) at 
different p values, L = AT = 16, A = 0, = 0, and 6g -O.St 
(figure), 6g = — 0.6t (inset). In both cases, macroscopic phase 
separation emerges ai p > pc- 

change with respect to the ^ = case is instead the 
size of the PS domains. Since the number of single 
atoms available for such a short-range antiferromagnetic 
arrangement diminishes with p, A should increase corre- 
spondingly, so that for p > Pc macroscopic phase segre- 



gation occurs also at ^ 7^ 0. This is indeed confirmed 
by numerical investigation of the charge structure fac- 
tor N{q) = e**^^ {{njfij-^r) — (^■)(^j+r)) at half chain 
{j = L/2). The results are reported in Fig. [3] at half 
filling. For p < Pc the value at which N{q) reaches its 
maximum moves smoothly from 27Tnd at p = to the low- 
est available q value (i.e., q = 27r/L) dit p = Pc, where it 
remains also at higher p values. In this case macroscopic 
phase segregation is achieved. The inset of Fig. [3] shows 
that the latter feature is maintained even in the crossover 
region g > —tad > though in this case the behavior of 
N{q) becomes closer to the standard Hubbard case for 

P < Pc 



IV. CONCLUSIONS 

In this paper, we provided exact analytical and numeri- 
cal evidence that a lattice Hamiltonian recently proposed 
[m, [1^ to mimic the physics of ultracold fermionic atoms 
on one-dimensional optical lattices close to a broad Fes- 
hbach resonance reproduces quite well the physics char- 
acteristic of these systems. In particular, aside to super- 
fluid and insulating phases, a phase-separated regime is 
achieved in case of strongly interacting atoms. In this 
case, macroscopic phase segregation appears in the den- 
sity profile above a critical polarization pc thanks to a 
peculiar mechanism. In fact, when PS is present in the 
system of balanced atoms, the coexisting phases form 
domains of microscopic dimension in order to implement 
antiferromagnetic short-range order. By increasing the 
imbalance of the population p, the size of the domains 
increases, to reach macroscopic phase segregation exactly 
at p = and above. The effect can be observed -for a 
large range of filling values around half filling- also on 
the number of pairs Ud- 

The model discussed here shows macroscopic phase 
segregation and Flude-Ferrell-Larkin-Ovchinnikov oscil- 
lations [17|] , both distinctive features in the investigation 
of systems of ultracold fermionic atoms. Therefore, our 
work gives evidence that Hubbard models with particle- 
assisted tunneling rates are more appropriate than the 
ordinary Hubbard model to describe the above physics 
when these systems are confined to anisotropic optical 
lattices. The present scenario should survive in higher 
dimension since its nature is to ascribe to the coopera- 
tive behavior of spin and charge degrees of freedom in 
the PS phase rather than to their separation. 



V. ACKNOWLEDGMENTS 

We thank F. Ortolani for the Lanczos code. This 
work was partially supported by the Grant No. PRIN 
2005021773 of Italian MIUR. 



5 



[1] M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and W. 

Ketterle, Science 311, 492 (2006). 
[2] M. W. Zwierlein, C. H. Schunck, A. Schirotzek, and W. 

Ketterle, Nature (London) 442, 54 (2006). 
[3] Y. I. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, 

Nature (London) 451, 689 (2008). 
[4] G. B. Partridge, W. H. Li, R. L Kamar, Y. A. Liao, and 

R. G. Hulet, Science 311, 503 (2006). 
[5] Y. Shin, M. W. Zwierlein, C. H. Schunck, A. Schirotzek, 

and W. Ketterle, Phys. Rev. Lett. 97, 030401 (2006). 
[6] S. Pilati and S. Giorgini, Phys. Rev. Lett. 100, 030401 

(2008). 

[7] D. Jaksch and P. Zoller, Ann. Phys. 315, 52 (2005). 

[8] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. 

De Sen, and U. Sen, Adv. Phys. 56, 243 (2007). 
[9] R. Jordens, N. Strohmaier, K. Gunther, H. Moritz, and 

T. Esslinger, Nature (London) 455, 204 (2008). 
[10] U. Schneider, L. Hackermuller, S. Wih, T. Best, L Bloch, 

T. A. Costi, R. W. Helmes, D. Rasch, and A. Rosch, 

Science, 322, 1520 (2008). 
[11] E. Timmermans, P. Tommasini, M. S. Hussein, and A. 

Kerman, Phys. Rep. 315, 199 (1999); K. Sengupta and N. 



Dupuis, Europhys. Lett. 70, 586 (2005); V. G. Rousseau 
and P. J. H. Denteneer, Phys. Rev. Lett. 102, 015301 
(2009). 

[12] J. K. Chin, D. E. Miller, Y. Liu, C. Stan, W. Setiawan, 

C. Sanner, K. Xu, and W. Ketterle, Nature (London) 

443, 961 (2006). 
[13] T. Stoferle, H. Moritz, K. Gunter, M. Kohl, and T. 

Esslinger, Phys. Rev. Lett. 96, 030401 (2006). 
[14] R. B. Diener and T.-L. Ho, Phys. Rev. Lett. 96, 010402 

(2006). 

[15] L.-M. Duan, Phys. Rev. Lett. 95, 243202 (2005). 
[16] L.-M. Duan, EPL 81, 20001 (2008). 

[17] B. Wang and L.-M. Duan, Phys. Rev. A 79, 043612 
(2009). 

[18] M. A. Simon and A. A. Aligia, Phys. Rev. B 48, 7471 
(1993). 

[19] See, for instance, M. Nakamura, T. Okano and K. Itoh, 
Phys. Rev. B 72, 115121 (2005) and references therein. 

[20] A. Montorsi, J. Stat. Mech. L09001(2008). 

[21] A. Anfossi, C. Degli Esposti Boschi, and A. Montorsi, 
Phys. Rev. B 79, 235117 (2009). 



